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Abstract. We present chemical abundances of a mild barium star HD202109 (£ Cyg) determined from the 
analysis of a spectrum obtained by using the 2-m telescope at the Peak Terskol Observatory and a high-resolution 
spectrometer with R = 80,000, signal to noise ratio >100. We also present the atmospheric parameters of the 
star determined using various methods including iron-line abundance analysis. For line identifications, we use 
whole-range synthetic spectra computed from Kurucz's database and the latest lists of spectral lines. Among the 
determined abundances of 51 elements, those of P, S, K, Cu, Zn, Ge, Rb, Sr, Nb, Mo, Ru, Rh, Pd, In, Sm, Gd, Tb, 
Dy, Er, Tm, Hf, Os, Ir, Pt, Tl, and Pb were not investigated previously. Assuming that the overabundance pattern 
of Ba stars is due to binary accretion, the observed abundance pattern of the neutron-capture process elements 
in HD202109 can be explained by combining the AGB star nucleosynthesis and the wind accretion scenario. 

Key words. Line: identification - Stars: abundances - Stars: atmospheres - Stars: evolution - Stars: chemically 
peculiar - Stars: binaries - Stars: individual: HD202109 



1. Introduction 

Barium stars were first denned by Bidclman & Keenan 
(1951). Burbidge & Burbidge (1957) obtained the first 
quantitative result about the chemical composition of a 
barium star HD46407 and reported overabundance of the 
heavy elements. Burbidge et al. (1957) interpreted the 
overabundance in terms of the s-process. Another impor- 
tant finding in understanding the physical process of the 
formation of barium stars was the fact that these stars 
are members of binary systems (McClure et al. 1980). 
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Wallerstein et al. (1997). revieved the development of the 
theory of element creation in stars. 

To explain the process of stellar evolution from abun- 
dance analysis, it is essential to construct detailed ob- 
served abundance patterns based on abundances of a 
large number of elements. Unfortunately, detailed abun- 
dance patterns are known only for few stars includ- 
ing the Sun, for which abundances of 73 chemical ele- 
ments are known (Grevesse & Sauval 1998). Some other 
well-determined stellar abundance patterns are those for 
GS22892-05 (based on 57 elements; Sneden et al. 2003), 
GS31082-001 (based on 38 elements with Z > 38; Aoki et 
al. 2003; Sneden et al. 2000), Procyon (based on 55 ele- 
ments; Yushchenko & Gopka 1996a,b), Przybylski's star 
(based on 54 elements; Cowley et al. 2000), and \ Lupi 
(based on 51 elements; Leckrone et al. 1999). 
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Theoretical abundance patterns of s-process elements 
in barium stars were published, for example, by Cowley & 
Downs (1980), Malaney (1987), Liang et al. (2000). The 
hypothesis that the overabundances of heavy elements in 
the atmospheres of barium stars are produced by accreting 
the ejecta from their AGB companion progenitor seems to 
be the most reliable explanation. These AGB companions 
have quickly evolved into white dwarfs and cannot be eas- 
ily detected. 

In this paper, we attempt to construct the detailed 
abundance pattern of one of the barium stars by deter- 
mining abundances of as many elements as possible, pay- 
ing special attention to s-process elements. We used high 
resolution observations in the blue spectral region, the lat- 
est atomic data, and the most updated spectrum synthesis 
method for better line identifications. 

The star selected for our investigation is HD202109 
(( Cyg, BS 8115, 64 Cyg), which is a prototype mild 
barium star. The star was first noticed as a barium star 
by Bidelman. The first quantitative abundance analysis 
for this star was done by Chromey et al. (1969). Keenan 
& MacNeil (1976) determined the star's spectral class as 
G8+III-IIIa Ba 0.6. Griffin & Keenan (1992) reclassified 
the spectral class as G8+IIIa Ba 0.6 and published the or- 
bit of the bright component of this spectroscopic binary. 
Griffin (1996) pointed out that the orbital period of the 
system (~ 18 years) is the longest among the systems with 
moderate enhancement of s-process elements. Pourbaix & 
Jorissen (2000) investigated the HIPPARCOS transit data 
of HD202109, but its long period did not permit them to 
find a reliable orbit. The white dwarf companion was di- 
rectly imaged in recent HST observations by Barstow et 
al. (2001). 

Previous studies on the chemical composition of the 
star include those of Sneden et al. (1981), Gratton (1985), 
Berdyugina (1993) (CNO elements), Zacs (1994), and 
Boyarchuk et al. (2001) (heavier elements). Cohen et al. 
(1999) and Garetta et al. (2001) used HD202109 as a com- 
parison star and published the abundances of several ele- 
ments. 

We select HD202109 as our target for the following 
reasons. First, since chemical abundances are known for 
some elements from previous studies, a comparison with 
the previous results allows us to derive more reliable abun- 
dance pattern. Second, the star is bright enough for us to 
acquire high-quality spectra using a mid-size telescope. 
Third, HD202109 is an ideal target to test the accretion 
hypothesis. Can we apply this hypothesis to barium stars 
with moderate enhancement of s-process elements, which 
is expected for barium stars with slow orbital motions? 
Fourth, since the star was observed in a wide range of 
wavelengths from X-ray to radio, we can derive a more 
concrete model of this binary system. 

2. Observations and data reduction 

A high resolution spectrum of HD202109 was obtained on 
September 19 in 2000 using a coude-echelle spectrometer 



(Musaev et al. 1999) mounted on the 2-m "Zeiss" telescope 
at the Peak Terskol Observatory located near Mt. Elbrus 
(Northern Caucasus, Russia) 3,124 m above sea level. 

We used one of the medium-resolution modes of the 
spectrograph with a resolving power of R = 80, 000. The 
observed wavelength range, AA3495 - 10000 AA, was cov- 
ered by 85 echelle orders. In the observed spectrum, there 
are gaps between the orders in the wavelength region of 
A > 3780 A and the width of each gap increases from 0.5 
A to 125 A as the wavelength increases. The CCD, made 
by Wright Instruments, has a 1,242 x 1, 152 format and 
each pixel has a physical size of 22 /im. The signal-to-noise 
ratio of the spectrum is not less than 100 over the whole 
wavelength range and reaches 300 or more in the blue and 
red parts of the spectrum. 

The first-stage data processing (background subtrac- 
tion, echelle vector extraction from the echelle- images, and 
wavelength calibration) was done using the latest ver- 
sion of PC-based DECH software (Galazutdinov 1992). 
For other processes including continuum placement, we 
use URAN software developed by one of us (Yushchenko 
1998). The location of the continuum was determined 
taking into account the calculated spectrum. More de- 
tails about this procedure were discussed for example by 
Walgren (1995). 

3. Atmospheric parameters 

With the processed data, we then determined the atmo- 
spheric parameters including the effective temperature, 
T e g, the surface gravity, log g, and the microturbulence, 

^micro- 

3.1. Effective temperature from colors 

We first determined T c g of HD202109 from color informa- 
tion. For this, we use two colors of V — K and b — y. 
Following the calibrations of Di Benedetto (1998) for 
V — K and Alonco et al. (1999) for b — y, we find that 
the temperatures determined from these two colors are 
T off = 4927 K and 4903 K, respectively. 

We compare our results of T e g with previous determi- 
nations based on different methods and data. Andrievsky 
et al. (2002) determined T e g of the star by using different 
photometric data and calibrations. They found five values 
of T cff =5130 K, 5070 K, 4900 K, 5120 K, and 5100 K and 
presented T e ff=5100 K as a representative value. Gray & 
Brown (2001) used the line-ratio method, and found the 
value T off =4987 K. Gray & Brown (2001) used HD202109 
as one of the calibrating stars to set their temperature 
scale. Upon our request, V. Kovtyukh kindly provided his 
determination based on the depth ratios of the iron lines. 
Using the latest version of the Kovtyukh & Gorlova (2000) 
method, he found T cB = 5044 ± 44 K. 

For the effective temperature one could use the tem- 
perature determined from color information, but we de- 
cided to choose the value determined from the iron-line 
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abundance analysis because it is known that the preci- 
sion of T e ff determined from spectroscopic determinations 
is significantly higher than that of the temperature deter- 
mined from color information (Boyarchuk 2001). 

3.2. Atmospheric parameters from iron line abundance 
analysis 

Usually atmospheric parameters from iron line analysis 
are determined by investigating the correlation(s) of the 
observed equivalent widths, log(WA/A), (and the excita- 
tion potential, E\ ow ) of iron lines with the iron abundances 
calculated based on the equivalent widths of the individ- 
ual lines by assuming a certain atmosphere model. In this 
process, it is customary to use fixed values of T e g and 
log g, under the assumption that these parameters can be 
constrained from other information (for example from col- 
ors), and leave v m icro as the only free parameter. Variation 
of this parameter makes it possible to avoid the influence 
of the uncertainties in effective temperature and surface 
gravity. 

Yushchenko et al. (1999) noticed that the scattering 
information of the iron abundances derived from differ- 
ent iron lines is useful in determining correct atmospheric 
parameters. The scattering for models with wrong param- 
eters is usually larger than for models with right parame- 
ters. 

In our analysis, we therefore determine the atmo- 
spheric parameters not only by leaving all of them as free 
parameters but also by using the additional information 
of scattering in the derived iron abundances. 

We selected 89 clean Fe I lines in the synthetic spec- 
trum of HD202109 and measured their equivalent widths 
in the observed spectrum by fitting their profiles with a 
gaussian. 

We used solar oscillator strengths (Gurtovenko & 
Kostik 1989). For iron lines where oscillator strengths 
were not available in Gurtovenko & Kostik (1989), we 
used solar values computed by the spectrum synthesis 
method. Details of the line selection and the synthetic 
spectrum calculations are described in Sect. 4. Once os- 
cillator strengths are set, the iron abundances are derived 
from the individual lines by using the WIDTH9 code of 
Kurucz (1995). 

HD202109 is a binary system. A full description of bi- 
nary systems requires the determination of two sets of 
atmospheric parameters for the individual components of 
the system and the flux ratio between them. We note, how- 
ever, that although HD202109 is a binary, the flux from 
the white dwarf can be neglected in our case. For determi- 
nation of atmospheric parameters of a binary where the 
fluxes of the component stars are of the same order, see 
Yushchenko et al. (1999). 

For atmosphere models Kurucz's (1995) data base was 
used. We made a 21 x 6 subgrid of atmosphere mod- 
els by subdividing the grid of Kurucz's data base. The 
ranges (and intervals) of the individual parameters are 
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Fig. 1. The correlation between the iron abundances de- 
termined from 89 Fe I lines in the spectrum of HD202109 
and the excitation energies, E\ ow , of the individual lines. 
Our final atmosphere parameters and abundances of other 
chemical elements were used. 



4750 K < T cff < 5250 K (AT cff = 25 K) for the effective 
temperature and 2.5 < log g < 3.0 (Alog g = 0.1) for 
the surface gravity. Calculations for all models were made 
with 41 values of the microturbulcnt velocity 1.0 km s _1 < 

^micro 

< 3.0km s 1 (Aw micro = 0.05 km s x ) 
Based on each set of atmosphere parameters, we com- 
pute iron abundances corresponding to the individual 
lines. For iron abundances computed for each set of at- 
mosphere parameters, we then calculate the scattering of 
the abundances and the coefficient (s) of the correlation(s) 
between the equivalent widths (and the excitation poten- 
tials) and the iron abundances. 

Among the tested models, the best atmosphere model 
was chosen to be the one providing zero (or very close to 
zero) correlation coefficients and minimal scattering in the 
derived iron abundances. We find that the best model has 
atmosphere parameters T c s — 5050 K, log g — 2.8, and 
I'micro = 1-45 km s . We note that our determination of 
T c ff is very close to the value of Kovtyukh determined by 
using the line-depth ratio method. We present the corre- 
lation between logiV(Fe) and \og{W\/\) in Fig.l and the 
correlation between logiV(Fe) and E\ OV! in Fig. 2 for the 
best atmosphere model. In Table 1, we compare our final 
atmospheric parameters with those of other determina- 
tions. 

Note that since our method uses scattering informa- 
tion of iron abundances, it enables us to determine the 
iron abundance. We found that the iron abundance of 
HD202109 is +0.01 ± 0.11 dex relative to the solar value. 



4. Abundance Analysis 

4.1. Synthetic Spectrum 

We calculated the abundances of all elements except iron 
by using the synthetic spectrum method. To construct ac- 
curate synthetic spectra one needs to find the parameters 
of line broadening. 

The magnetic field of HD202109 was investigated by 
Tarasova (2002). She found values of -5.4±0.2, -0.4±2.0 
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Fig. 2. The correlation between the iron abundances de- 
termined from 89 Fe I lines in the spectrum of HD202109 
and the equivalent widths, \og{W\/X), of the individual 
lines. 



Table 1. Atmospheric parameters of HD202109 (C Cyg) 



reference 


T c s 


log g 


^micro 




(K) 




(km s" 1 ) 


Chromey (1969) 


5143 






Pilachowski (1977) 


4893 


2.8 




Sneden et al. (1981) 


4870 


2.5 


2.0 




5000 


2.8 


2.2 


Gratton et al.(1982) 


4941 


2.0 




Gratton (1985) 


4950 


2.0 


2.4 


Fernandez- Villacanas 


4900 


2.0 




et al. (1990) 








McWilliam (1990) 


4990 


2.87 




Berdyugina (1993) 


5000 


2.7 


2.0 


Berdyugina & Savanov (1994) 


5000 


2.8 




Zacs (1994) 


5050 


2.8 


3.5 


Cohen et al. (1999) 


4950 


2.7 


1.6 


Boyarchuk et al. (2001) 


4977 


2.52 


1.4 


Gray & Brown (2001) 


4987 






Andrievsky et al. (2002) 


5100 


2.5 


1.5 


This paper 


5050 


2.8 


1.45 



and 0.2±1.8 gauss. A field of this strength cannot influence 
our final final abundances significantly. The rotation veloc- 
ity was determined by Gray (1989) v mt — 3.4 ± 0.5km s _1 
and Medeiros & Mayor (1999) v rot ^ 1.0 km s -1 . 

We assume that lines are broadened mainly by the 
macroturbulence. We estimate the macroturbulent veloc- 
ity by analyzing the profiles of iron lines with accurate os- 
cillator strengths. We find that t> ma cro = 3.2 — 3.5 km s~ . 
This value includes all possible line broadening mecha- 
nisms. We used a Gaussian model of macroturbulence 
(Gray 1976). 

With the input of all these parameters, the syn- 
thetic spectra are constructed using the SYNTHE code 
of Kurucz (1995). We used various lists of spectral lines. 
These include all atomic and molecular lines of Kurucz's 
database (1995), Morton's (2000) lines, the DREAM 



database (Biemont et al. 2002) lines for lanthanides and 
actinides, part of the lines from VALD database (Piskunov 
et al. 1995), and lines in other lists. We took into account 
the hyperfine-structure lines and isotopic splitting of the 
lines for Mn, Cu, and Eu. Split-line data are taken from 
the Kurucz database (1995). 

4.2. Line identifications 

To better identify heavy element lines, the synthetic spec- 
trum used for line identifications is produced by increasing 
the abundances of uninvestigated heavy elements by an 
amount of +0.5 dex compared to the solar system values. 

For identification of lines of any element we selected all 
not-strongly-blended lines of this element in the results 
of synthetic spectrum calculations and inspected these 
lines in the observed spectrum to select the best lines for 
abundance analysis. The observed and synthetic spectra of 
HD202109 and of the Sun were displayed simultaneously 
on the computer screen. This permits to avoid identifica- 
tion errors and to increase the number of selected lines. 
The URAN code of Yushchenko (1998) was used. 

In Figs. 3 and 4, we present a part of the observed spec- 
trum of HD202109 and its approximation by the synthetic 
one. 

4.3. Abundance Determinations 

Once the individual lines are identified, we determine 
the element abundances. For line identifications and 
abundance determinations we used the URAN code of 
Yushchenko (1998). In this code abundances are com- 
puted in semiautomatic mode. To approximate the ob- 
served spectrum by a synthetic one we changed the os- 
cillator strengths of all lines in the investigated region, 
except for those lines lying within 0.03 A of the investi- 
gated line. The SYNTHE code of Kurucz (1995) was used 
for calculations of the synthetic spectra. 

More descriptions of this type of codes for semiau- 
tomatic and automatic abundance analysis are found in 
Cowley (1995), Tsymbal & Cowley (2000) (MERSEN 
code), Valenty & Piskunov (1996) (SME), Erspamer & 
North (2002), and Bruntt et al. (2002). Although these 
codes use different algorithms, they allow one to signifi- 
cantly increase the list of investigated lines and minimize 
manual work. 

As initial abundances for the computation of synthetic 
spectrum, we use the values published in earlier investi- 
gations. This permits us to include moderately blended 
lines for abundance analysis and thus maximize the num- 
ber of elements with determined abundances. In addition, 
a good guess of the initial abundances can minimize the 
computation time by reducing the number of iterations in 
the fitting process. 

To compare the abundances of individual elements of 
HD202109 with those of the Sun, we also calculate the 
solar abundances for the investigated lines. For this calcu- 
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_N %_ 

11 _915 
12_884 
13_943 
14_872 
15_562 
1B_953 
17 1000 
18_549 
19_926 

20 _787 

21 _782 

22 _587 
23_619 

24 _564 

25 _726 

26 _561 

27 _546 

28 _885 

29 _808 

30 _988 . 

31 _845 

32 _635 . 

33 _627 

34 _726 . 

35 _871 

36 1000 

37 1000 

38 548 



_R Wave Id 

471 32.328 Nd II = 
825 32.490 C2 
477 32.669 Fe II 
900 33.158 MgH 
948 33.239 CN 
602 33.466 Co I 
111 33.681 Fe I 
721 33.823 C2 
816 34.200 MgH 
800 34.313 C2 
808 34.482 C2 
700 34.558 MgH 
770 34.666 C2 
803 35.072 MgH 
787 35.140 Pr II # 
711 35.570 C2 
708 35.578 C2 
823 35.684 C2 
877 35.945 Cr I 
489 36.091 Fe I 
820 36.270 C2 
772 36.437 C2 
884 36.549 Ru I = 
812 36.658 C2 
593 36.802 Fe II 
.105 37.070 Nl I 
152 37.395 Fe I 
698 37.584 C2 




Fig. 3. Part of the spectrum of HD202109. The squares 
are the observed spectrum. The solid line is the synthetic 
one. The positions of the strong and faint lines in the syn- 
thetic spectrum are shown by long and short dashes at the 
bottom of the figure. Part of the strong lines are marked 
by numbers and identification. A table with line data is in 
the left of the figure. The first column is the line number. 
The second column is the portion of the line in the total 
line absorption coefficient at the wavelength of the center 
of the line in the synthetic spectrum. For a clean line the 
value in this column must be 1000. The third column is 
the value of synthetic spectrum at the center of the line. 
The continuum value is 1000. Only strong lines are listed. 
The values of the synthetic spectrum in the table are not 
smoothed by instrumental and macroturbulence profiles. 
In the last two columns we give the last digits of the wave- 
lengths and the identifications. Lines of r-, s-process ele- 
ments are marked by an equal sign. The Pr II A5135.140 
line is marked by a # sign in the table and by a vertical 
dashed line in the spectrum. This and the next figures are 
the PrintScreen output of URAN software (Yushchenko 
1998). 

lation, we use the Liege Solar atlas (Delbouille et al. 1974) 
and the Holweger & Muller (1974) atmosphere model. The 
adopted values of the microturbulent and macroturbu- 
lent velocities are 1.0 km s _1 (Gopka & Yushchenko 1995) 
and 1.8 km s~ , respectively. The continuum in the Liege 
Solar atlas is corrected in accordance with Arderberg & 
Virdeforce (1979) and Rutten & van der Zalm (1984). We 
use the same SYNTHE and URAN codes for synthetic 
spectrum construction, line identifications, and abundance 
determinations . 

We present a partial list of the determined abundances 
of elements derived from the individual lines of HD202109 
in Table 2 (for iron lines) and 3 (for lines of other ele- 
ments). Table 2 contains the element code (26.00 for Fel), 
wavelength, equivalent width (Eq.W.), oscillator strength 
(log gf), excitation potential, and abundance (expressed 
on the scale of logAf(H) = 12). 



_N %_ 

1 1 _054 
12_078 
13_999 
14_237 
15_356 
16_181 
17_180 
18_102 
19_828 
20_216 

21 _849 

22 _267 

23 _444 

24 _496 
25_147 

26 _075 

27 _846 
28_152 

29 _051 

30 _295 

31 _086 

32 _954 

33 _069 

34 _056 
35_189 

36 _347 

37 _058 

38 450 



_R Wave Id 

919 59.179 MgH 
734 59.200 CH 
223 59.272 Ti I 
975 59.388 MgH 
979 59.428 CN 

979 59.447 MgH 

980 59.463 Cr I 
980 59.468 CN 
957 59.557 CH 
800 59.627 Ti I 
747 59.653 Erll# 
868 59.714 CH 
877 59.726 Cr I 
931 59.772 CH 
967 59.838 MgH 
965 59.840 CN 
787 59.912 Cr I 
791 59.926 Ce II = 
910 59.974 MgH 
884 60.000 CH 
588 60.048 Dy II = 
542 60.068 Fe I 
926 60.153 CN 
926 60.156 CN 

920 60.169 Fe II 
913 60.182 Mo I = 
912 60.189 MgH 
919 60.205 Ni I 
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Fig. 4. Observed (squares) and calculated spectra of 
HD202109. The solid line is the spectrum, calculated with 
our final abundances. Dashed lines are the spectra, cal- 
culated with erbium abundances changed by ±0.2 dex 
from the best value. The central line of the plot ( Er II 
A4759.653 ) is marked by a # sign in the table and by a 
vertical dashed line in the spectrum. 

Table 2. Fe I and Fe II lines in the spectrum of HD202109 
(example, full table is available in electronic form) 



element 


A 


Eq.W. 


log gf 


-Blow 


logiV 


code 


(A) 


(mA) 




(eV) 




26.01 


4620.52 


85 


-3.57 


2.828 


7.827 


26.01 


5197.58 


114 


-2.48 


3.230 


7.693 


26.01 


5234.63 


115 


-2.42 


3.221 


7.635 


26.01 


6369.46 


43 


-4.37 


2.891 


7.707 


26.01 


6432.68 


71 


-3.83 


2.891 


7.728 


26.01 


6456.38 


85 


-2.32 


3.903 


7.594 



Table 3 contains the element code, wavelength, oscil- 
lator strength, excitation potential, abundances for £ Cyg 
and for the Sun and the difference between these abun- 
dances (A). For lines without counterparts in the solar 
spectrum Grevesse & Sauval (1998) solar system abun- 
dances were used. The full table in electronic form can 
be obtained at the Web sites: "users . Odessa. net/~yua" 
and "yushchenko.netfirms.com" 

In Tables 4 and 5, we present the final list of the abun- 
dances of 51 elements, which are computed by averaging 
the abundances derived from all lines of the individual el- 
ements. Table 4 contains CNO data obtained by Sneden 
et al.(1981), Gratton (1985), Berdyugina (1993), Cohen et 
el. (1999), and our result for oxygen. 

The relative abundances of the other elements (listed 
in Table 5) are compared with the determinations of Zacs 
(1994) and Boyarchuk et al. (2001). Our results (except 
iron) were obtained with the spectrum synthesis method. 
The Zacs (1994) and Boyarchuk et al.(2001) data were 
calculated with the model atmospheres method. Elements 
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Table 3. Abundances of chemical elements calculated 
from individual lines in the spectrum of £ Cyg (HD202109) 
and in the solar spectrum (example, full table is available 
in electronic form) 



element 


A 


i°g gf 


-C'low 


1 ^ 7V 7" 

lOgiV 




A 
A 


code 


(A) 




(pV) 

l eV ) 


r C'vtr 


O {111 




72.01 


4093.155 


-1.090 


.452 


1.01 


.56 


.44 


76.00 


3528.598 


-1.700 


.000 


1.46 


1.41 


.05 


76.00 


3501.163 


-.510 


1.841 


1.96 


1.41 


.55 


77.00 


3515.947 


-1.580 


.881 


< 1.75 






77.00 


3594.388 


-1.550 


.881 


< 1.85 






78.00 


3638.789 


-1.340 


1.254 


1.8 






81.00 


3519.210 


.141 


.966 


< 1.4 






82.00 


3683.462 


-.555 


.969 


< 2.1 






82.00 


3639.567 


-.715 


.969 


< 2.1 







with lines having no counterparts in the solar spectrum are 
marked by the asterisk. For these elements abundances are 
given with respect to the solar values (Grevesse & Sauval 
1998) of the corresponding elements. Also listed are the 
numbers of lines, (n), from which the mean abundances 
of the individual elements are computed. The uncertainty 
of each element's abundance is estimated by computing 
the standard deviation of the abundances derived from 
the individual lines of the clement. The uncertainties are, 
therefore, presented only for elements whose abundances 
are based on more than two lines. 

We investigated the abundances of 47 elements. The 
abundances of Li, C, N, and Ba are taken from the papers 
mentioned above. The total abundance sample consists of 
51 elements, s-process elements show overabundances in 
the atmosphere of this star. 

The abundances of P, S, K, Cu, Zn, Ge, Rb, Sr, Nb, 
Mo, Ru, Rh, Pd, In, Sm, Gd, Tb, Dy, Er, Tm, Hf, Os, Ir, 
Pt, Tl, and Pb have been determined for the first time in 
this paper. 

4.4. Uncertainties in abundances 

The abundances of elements are subject to uncertainties 
caused by various sources. The uncertainties in the deter- 
mined abundances caused by the uncertainty of the atmo- 
sphere model can be estimated by investigating how the 
iron abundance varies depending on the adopted values 
of atmospheric parameters. We find that the variations 
of the iron abundance caused by the small deviations of 
AT off = 100 K, Aw micro = 0.1 km s"\ and Alog g = 0.2 
are 0.06 dex, 0.10 dex, and 0.01 dex, respectively. 

Of course, the abundances of different species will have 
different dependencies on the perturbations in the atmo- 
sphere parameters. For example, we find that the car- 
bon abundance estimated based on the atmosphere model 
adopted by Cohen et al. (1999) differs from our determi- 
nation by an amount of 0.23 dex. We note, however, that 
our adopted atmospheric parameters of HD202109 agree 
well with other determinations and thus we believe that 



Table 4. CNO abundances in the atmosphere of 
HD202109 with respect to their abundances in the solar 
atmosphere 





Sneden 


Gratton 


Berdyu- 


Cohen 


This 




et al. 




gina 


et al. 


work n 




(1981) 


(1985) 


(1993) 


(1999) 




c 


-0.18 


-0.10 


-0.19 






N 


+0.58 


+0.61 


+0.23 






O 


+0.04 


-0.34 


+0.02 


+0.35 


-0.22: 1 



the uncertainty in the determined abundances due to the 
adopted atmosphere model is unlikely to be large. 

If we adopt one of the parameters with some uncer- 
tainty, we try to select the other parameters to minimize 
correlation coefficients and scattering of the results. For 
example, the metallicity derived in our investigation is 
+0.01 dex, and if we adopt the parameters of Boyarchuk 
et al. (2001) it will be changed by 0.04 dex only. 

5. Observed abundance pattern 

5.1. CNO elements 

From Table 4, one finds that the abundances of these ele- 
ments show considerable differences from one determina- 
tion to another. In our analysis, we determine the abun- 
dance of oxygen only. We investigate the cause of the dif- 
ferences, focusing mainly on the difference between the 
oxygen abundances determined by us and by Cohen et al. 
(1999), for which the difference is greatest. We find several 
possible causes of differences. 

First, the analyses are based on different oxygen lines. 
Our oxygen abundance is derived from line A6300.304. 
On the other hand, the oxygen abundance determined by 
Cohen et al. (1999) was based on triplet AA7771-7775. 
Since this triplet is located in a gap in our spectra, we 
could not use it for the abundance analysis. 

Second, the lines used for both analyses are affected 
by blending. The A6300.304 line used for our analysis is 
contaminated by the nickel line. Among the three lines 
used by Cohen et al., only the A7771.944 line is free from 
blending and the other two lines are contaminated by 
CN lines. The difference in oxygen abundances, calculated 
with the Cohen et al. (1999) value of the equivalent width 
of A7771.944 and calculated using spectrum synthesis for 
A6300.304 line is only 0.09 dex. 

Third, the difference between the adopted atmosphere 
models also has some influence on the oxygen abundances. 
We calculated the oxygen abundance from A6300.304 line 
based on the atmosphere model adopted by Cohen et al. 
(1999). This abundance is higher than our original oxygen 
abundance by an amount of 0.23 dex. 

CNO abundances are important for the abundance 
determinations of other elements. We calculated several 
synthetic spectra for the whole observed region with our 
atmosphere model parameters and different CNO abun- 
dances. We chose to use the values of Gratton (1985). 
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Table 5. The abundances of chemical elements (except 
abundances in the solar atmosphere. 



[N/Nh] (n) 



element 


ZelCS 

(1994) 




Boyarch- 
uk et al. 
(2001) 


This work 


3 Li I 


—0.14 


(11 










8 O I 








—0.22: 




(11 
y ± j 


11 Na I 


—0.35 


(21 
y^j 


+0.19 


+0.24+ 


.08 


(41 
y^j 


12 Mg I 


—0.51 


(21 




+0.22+ 


22 


(61 
y u ) 


13 Al I 






+0.16 


+0.12+. 


12 


(81 
y°t 


14 Si I 


+0.14±.13 


(3) 
y^j 


+0.09 


—0.05+ 


.13 


(521 


Si II 








+0.26 




(2) 

y^> 


15 P I 








+0.10 




(11 
y ± j 


16 S I 








+0.00+ 


.12 


(31 

y°j 


19 K I 








—0.17 




(2) 
y^j 


20 Ca I 


+0.08±.19 


(4) 
y^j 


-0.03 


+0.02+ 


.09 


(51 
y°> 


21 Sc I 


+0 04± 30 


(41 
y^i 


—0.02 


+0.06+ 


.11 


(4) 
y^/ 


Sc II 








+0.08+ 


.15 


dol 


22 Ti I 


— 0.17±.22 


(21) 


—0.11 


—0.20+ 


.17 


(701 


Ti II 








—0.03+ 


.13 


(301 


23 V I 


— 0.13±.18 


(161 


—0.04 


—0.03+ 


.14 


(41) 


V II 








+0.17 




(1) 


24 Cr I 


—0 06± 18 

V/ . V 7 V ' — 1 — > 1 v 1 


fill 


—0.08 


-0.17+ 


.13 


(65) 


Cr II 








+0.10+ 


.10 


(15) 


25 Mn I 


—0 30± 13 


(51 
y°j 




-0.25+ 


.17 


(21) 


26 Fe I 


+0.12±.23 


(51) 


-0.03 


+0.01+ 


.11 


(89) 


Fe II 








+0.06+ 


.07 


(6) 


27 Co I 


-0.22±.12 


(6) 


-0.13 


+0.07+ 


.11 


(28) 


28 Ni I 


-0.05+.25 


(10) 


-0.09 


-0.09+ 


.17 


(74) 


29 Cu I 








-0.01 




(1) 


30 Zn I 








-0.08+ 


.09 


(4) 


32 Ge I 








+0.08 




(1) 


37 Rb I* 








-0.07 




(1) 


38 Sr I 








+0.26 




(2) 



Cratton (1985) abundances fit the observed spectra best 
among the determinations listed in Table 4. 

5.2. Na to K 

Our results for these elements are close to the Boyarchuk 
et al. (2001) data for common elements. On the other 
hand, our values show some differences from the data of 
Andrievsky et al. (2002), which were released just after 
we completed the major part of our abundance deter- 
minations. For example, the NLTE value of the sodium 
abundance determined by Andrievsky et al. (2000) is 
+0.35 dex, compared to our value of +0.24 dex. We sus- 
pect that the differences may be caused by their adoption 
of an NLTE atmosphere model. 

Our Mg abundance is in a good agreement with that 
of Cohen et al. (1999). Although the abundances of Si 
determined by Boyarchuk et al. (2001) and Cohen et al. 
(1999) are slightly higher than our value, all these values 
including ours are consistent with the solar abundance 
within the errors. Our results for P, S, and K indicate 
that the abundances of these elements are close to the 
solar values. 
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CNO) in the atmosphere of HD202109 with respect to the 



[N/Nu] (n) 



element Zacs Boyarch- This work 

(1994) uk et al. 
(2001) 



39 


Y I 


+0.37+.24 


(31 


+0.30 


+ 0.15+.15 


(3) 
y j j 




Y II 








+ 48+ 16 


(22) 


40 


Zr I 


— 0.08+.20 


(51 
y°j 




+ 0.21+.06 


(ID) 




Zr II 








+ 0.52+.18 


(12) 


41 


Nb II* 








+ 0.43 


(1) 


42 


Mo I 








+ 0.20 


(2) 
y^j 


44 


Ru I* 








4- 32+ 08 


(3) 
y°j 


45 


Rh I* 








< +0.2 


(2) 

y^j 


46 


Pd I* 








+ 0.32: 


(l) 


49 


In I 








+ 0.11: 


(1) 


56 


Ba II 


+0.41+.13 


(31 
y°i 


+0.54 






57 


La II 


+0.38+.15 


(3) 
y^j 


+0.45 


+ 0.51+.20 


(12) 


58 


Ce II 


+0.55 


(1) 


+0.33 


-1- n 37+ is 


(45) 


59 


Pr II 


+0.32+.13 


(3) 


+0.43 


+ 0.19+.19 


(6) 


60 


Nd II 






+0.23 


+ 0.42+.17 


(701 


62 


Sm II 








-1- n 34+ 18 


(151 


63 


Eu II 


-0.05 


(2) 


+0.22 


+ 0.32+.12 


(3) 


64 


Gb II 








+ 0.27+.19 


(4) 


65 


Tb II* 








+ 0.1 


(1) 


66 


Dy II 








+ 0.33+.18 


(4) 


68 


Er II 








+ 0.35 


(1) 


69 


Tm II 








< +0.2 


(1) 


72 


Hf II* 








+ 0.45: 


(1) 


76 


Os I* 








+ 0.30: 


(2) 


77 


Ir I* 








< +0.4 


(2) 


78 


Pt I* 








+ 0.0 


(1) 


81 


TI I* 








< +0.5 


(1) 


82 


Pb I* 








< +0.2 


(2) 



5.3. Ca to Ni 

These elements produce the majority of lines in the spec- 
tra of normal stars. We find that our iron abundance, 
+0.01 dex, agrees well with most previous determina- 
tions, ranging from -0.17 to +0.10 dex (Cayrel de Strobel 
et. al. 1977). Only Andrievsky et al. (2002) result is 
slightly higher: +0.15 dex. We suspect that the higher 
value of Andrievsky et al. (2002) is caused by the use of 
narrow-wavelength-range spectra in their abundance anal- 
ysis. Because of this their [Na/Fe] for HD202109 can be 
higher. 

We also find that the abundances derived from neutral 
and ionized lines result in similar values for Fe and Sc, 
while the values of Ti, V, and Cr show differences of ~ 
0.2 dex. 

The abundance of Mn, which is estimated by using hfs 
data from the Kurucz (1995) data base, indicates that this 
element is slightly underabundant. 

For elements in this group, we find that the determined 
abundances show a large fluctuation compared to the solar 
values. These fluctuations may be an intrinsic characteris- 
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Fig. 5. The abundances of chemical elements and ions in 
the atmosphere of HD202109 with respect to their abun- 
dances in the solar atmosphere. Our data and data of 
Boyarchuk et al. (2001), Zacs (1994), and Gratton (1985) 
are shown for comparison. 

tic of the abundance patterns of HD202109, or caused by 
inadequate sophistication of the atmospheric model. 

5.4. r-, s-process elements 

For Y and Zr, we find inconsistencies between the abun- 
dances derived from the lines of neutral and ionized 
species. The mean result for two ions of Y is in perfect 
agreement with Boyarchuk et al. (2001) and Zacs (1994). 

For line identifications of lanthanides, we use both the 
VALD line list (Piskunov et al. 1995) and the DREAM 
(Biemont et al. 2002) data base. We expected that by 
using the DREAM data base we could identify more lan- 
thanide lines based on the previous experience in the spec- 
tral analysis of Przybylski's star (Yushchenko et al. 2002). 
For HD202109, however, we could identify only a few more 
lines. We find that our results are consistent with the data 
of Boyarchuk et al. (2001) and Zacs (1994) for common el- 
ements. 

From our investigation of the spectrum of HD202109, 
we are able to determine abundances of r-, and s-process 
elements whose abundances were not previously known. 
These elements include Cu, Zn, Ge, Rb, Nb, Mo, Ru, Rh, 
Pd, In, Tb, Er, Tm, Hf, Os, Ir, Pt, Tl, and Pb. For most 
of these elements, we calculated the abundances based on 
one or two lines and for some we could only set upper lim- 
its. However, even this information is useful in constrain- 
ing the parameters of the wind accretion model, which is 
described in the next section. 

6. Predicted heavy element abundances 

As presented and shown in Table 5 and Fig. 5, respectively, 
neutron-capture (n-capture) process elements (hereafter 



25 35 45 55 65 75 85 

Z 

Fig. 6. The comparison of the observed composition of 
heavy elements in the atmosphere of HD202109 with the 
predicted abundances of these elements. Solid line - calcu- 
lation with a=1.5, dashed line - a— 1.8. Circles and trian- 
gles - observational data. For all elements except barium 
we plotted our results from Table 5, the Boyarchuk et al. 
(2001) value was plotted for barium. Big circles - more 
than 1 line were measured for this element, midium-sized 
circles - one line with respect to Sun, small circles - one 
line without direct comparison with the solar spectrum. 
Triangles - upper limit. For Y and Zr the weighted mean 
of the abundances from the lines of first and second spec- 
tra are plotted. 



heavy elements) of HD202109 are overabundant. Two 
peaks at Z^39-40 and Z^56-57 are obviously due to the 
neutron magic numbers of 50 and 82 nuclei in the path 
of s-process nucleosynthesis occurring in the interiors of 
AGB stars. 

What causes the heavy-element overabundances in 
barium stars? It is generally believed that the overabun- 
dances were caused by binary accretion, where Ba stars ac- 
creted the ejected material from their companions, the for- 
mer AGB stars and the current white dwarfs, which syn- 
thesized these heavy elements by themselves and ejected 
the elements into the interstellar medium through the stel- 
lar wind (Liang et al. 2000 and references therein). 

Liang et al. ( 2000 ) calculated the s-process nucleosyn- 
thesis of AGB star with mass 3M Q and solar metallic- 
ity. At the same time, they set up an angular momentum 
conservation model of wind accretion in binary systems. 
Combining the AGB star nucleosynthesis and the wind ac- 
cretion model, the theoretical heavy-element abundances 
of barium stars were calculated, and the observed abun- 
dances of some sample barium stars were explained suc- 
cessfully. 
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In this section, we check whether the binary accretion 
model can also explain the heavy element overabundance 
of HD202109 by comparing the abundance pattern deter- 
mined from our spectral analysis with the theoretical one 
calculated based on the model of Liang et al. (2000). 

Following Liang et al. (2000), we compute the theo- 
retical abundances of heavy elements in two steps. In the 
first step, we calculate the overabundances of the intrin- 
sic AGB star at each ejection by adopting the theory of 
s-process nucleosynthesis and the latest TP- AGB model 
(Straniero et al. 119951 Straniero et al. 119971 Gallino et al. 
1998] Busso et al. H999f) . In the second step, the overabun- 
dances of heavy elements in the atmosphere of barium star 
are calculated by accreting the ejected matter predicted 
from the model of wind accretion and mixing them on suc- 
cessive occasions. The more details about the scenarios of 
s-process nucleosynthesis and the orbital evolution of the 
binary system can be found in Liang et al. ( 2000 ) and Liu 
et al. iPTUUjl . 

The second step needs the main-sequence masses of 
the stars. Thus, we estimated the mass of HD202109 by 
using the stellar evolution tracks given by Girardi et al. 
(2000). The derived mass from the Mboi — logT c ff diagram 
is ~ 3.05 M Q , which reveals that HD202109 is a mild 
barium star (Jorissen et al. 1998; Liang et al. 2003). 

Therefore, the adopted wind accretion model is: 
3.0 Mq and 2.5 Mq for the main-sequence masses of the 
intrinsic AGB star (the current white dwarf companion) 
and the barium star, respectively; v c - } — 15 km s _1 for the 
wind velocity; 0.15 times the Bondi-Hoyle accretion rate 
for the actual accretion rate (Liang et al. 2000). 

The obtained orbital eccentricity and period of the bar- 
ium star system match the values e = 0.22 and P — 6489 
days of HD202109, respectively, which were determined 
observationally by Griffin & Keenan (1992). 

Fig. 6 gives the theoretical abundance pattern of 
HD202109 (the lines) and the pattern determined from 
our spectral analysis (the points). The two theoretical 
abundance pattern curves are the results of the cases with 
neutron exposures a — 1.5 (solid curve) and 1.8 (dashed 
curve), "a" refers to the neutron exposure occurring in 
its AGB star companion, and represents the times of the 
corresponding exposures in the 13 C profile suggested by 
Gallino et al. (1998). The more details can be found in 
Liang et al. (2000). 

From Fig. 6, one finds that the calculated abundance 
pattern is consistent with the detailed pattern of the ob- 
served abundances, which means that the heavy element 
overabundances of HD202109 may result from accreting 
the ejecta of its AGB star companion. The neutron ex- 
posure of the s-process nucleosynthesis occurring in the 
interior of the AGB companion corresponds to a — 1.5- 
1.8. 



7. Conclusion 

In this paper we make a detailed analysis of abundances of 
s-process elements in the atmosphere of the mild barium 
star HD202109. 

The atmosphere parameters were found from the care- 
ful analysis of iron abundances, calculated from individual 
iron lines. Taking into account the variation of scattering 
of the mean iron abundance for atmosphere models with 
different parameters, we were able to develop a method to 
determine the effective temperature, surface gravity, mi- 
croturbulent velocity and metallicity from the analysis of 
the iron lines. 

We used differential spectrum synthesis to find the ele- 
ment abundances, taking the solar spectrum as a compar- 
ison one. The result is the abundances of 47 elements in 
the atmosphere of HD202109. We found the abundances 
of several elements near the last peak of the abundances 
of s-process elements in barium stars - those near atomic 
number Z=80. The abundances of Li, C, N, Ba were known 
from previous investigations of this star. The total abun- 
dance sample consists of 51 elements. It is the most de- 
tailed abundance pattern for barium stars now available. 

We calculated the theoretical abundances of heavy ele- 
ments of barium stars using AGB star nucleosynthesis and 
the wind accretion model in barium binary systems. The 
observed abundances of heavy elements of HD202109 are 
consistent with the predicted abundances. It was shown 
that the barium star HD202109 can be formed through a 
wind accretion scenario. The corresponding neutron expo- 
sure in the s-process nucleosynthesis in the interior of its 
AGB star companion is a = 1.5 to a = 1.8. 

Barium stars with orbital period P >1600 days can be 
formed through wind accretion (Liang et al. 2000 Liang et 
al. 12)03). Jorissen et al. (1998) suggested that the corre- 
sponding period is 1500 days. Possibly, barium stars with 
lower orbital period form through other scenarios: dynam- 
ically stable late case C mass transfer or common envelope 
ejection. HD202109 (P=6489 days) provides strong sup- 
port for this suggestion by the consistency between the 
observed abundances of so many heavy elements and the 
corresponding predicted abundances from our wind accre- 
tion model. 
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